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FEYNMAN-KAC FORMULA FOR LEVY PROCESSES 
WITH DISCONTINUOUS KILLING RATE 

KATHRIN GLAU 


Abstract. The challenge to fruitfully merge state-of-the-art tech¬ 
niques from mathematical finance and numerical analysis has in¬ 
spired researchers to develop fast deterministic option pricing meth¬ 
ods. As a result, highly efficient algorithms to compute option 
prices in Levy models by solving partial integro differential equa¬ 
tions have been developed. In order to provide a solid mathe¬ 
matical foundation for these methods, we derive a Feynman-Kac 
representation of variational solutions to partial integro differen¬ 
tial equations that characterize conditional expectations of func¬ 
tionals of killed time-inhomogeneous Levy processes. We allow for 
a wide range of underlying stochastic processes, comprising pro¬ 
cesses with Brownian part as well as a broad class of pure jump 
processes such as generalized hyperbolic, multivariate normal in¬ 
verse Gaussian, tempered stable, and a-semi stable Levy processes. 

By virtue of our mild regularity assumptions as to the killing rate 
and the initial condition of the partial integro differential equation, 
our results provide a rigorous basis for numerous applications, in 
financial mathematics, probability theory and physics. We reen¬ 
counter the original ideas of Feynman and Kac, but now revealing 
the normal inverse Gaussian process in its role connecting the rel¬ 
ativistic Schrodinger equation to stochastic processes. In Regard 
to finance we suggest a flexible class of employee options. We im¬ 
plement a Galerkin scheme to solve the attendant pricing equation 
numerically and illustrate the effect of a killing rate. 

Time-inhomogeneous Levy process, killing rate, Feynman-Kac representation, weak 
solution, variational solution, parabolic evolution equation, partial integro differ¬ 
ential equation, pseudo differential equation, nonlocal operator, fractional Laplace 
operator, Sobolev-Slobodeckii spaces, option pricing, Laplace transform of occupa¬ 
tion time, relativistic Schrodinger equation, employee option, Galerkin method 
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1. Introduction 

Feynman-Kac formulas play a distinguished role in probability theory and func¬ 
tional analysis. Ever since their birth in 1949, Feynman-Kac-type formulas have 
been a constant source of fascinating insights in a wide range of disciplines. They 
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originate in the description of particle diffusion by connecting Schrodinger’s equa¬ 
tion and the heat equation to the Brownian motion, see (33)- A type of Feynman- 
Kac formula also figures at the beginning of modern mathematical finance: In their 
seminal article of 1973, Black and Scholes derived their Nobel Prize-winning op¬ 
tion pricing formula by expressing the price as a solution to a partial differential 
equation, thereby rediscovering Feynman and Kac’s deep link. 

The fundamental contribution of Feynman-Kac formulas is to link stochastic 
processes to solutions of deterministic partial differential equations. Thus they 
also establish a connection between probability theory and numerical analysis, two 
disciplines that have evolved largely separately. Although both enjoy great suc¬ 
cess, transfer between them has remained only incidental. This may very well be 
the reason for applications of Feynman-Kac still appearing so surprisingly fresh. 
In computational finance, they enable the development of option pricing methods 
by solving deterministic evolution equations. These have proven to be highly ef¬ 
ficient, particularly when compared to Monte Carlo simulation. Thus, like other 
deterministic methods, they come into play whenever efficiency is essential and 
the complexity of the pricing problem is not too high. This is the case for recur¬ 
ring tasks, such as calibration and real-time pricing, and over the last few decades 
has given rise to extensive research in computing option prices by solving partial 
differential equations. The challenge to extend these methods to price options in 
advanced jump models has furthermore inspired researchers in recent years to de¬ 
velop highly efficient and widely applicable algorithms, see for instance (|14r ). Hilber, 
Reich, Schwab and Winter (2009), Hilber, Reichmann, Schwab and Winter (2013), 
Salmi, Toivanen and Sydow (2014) and (f28i). 

In this article we derive a Feynman-Kac-type formula so as to provide a solid 
mathematical basis for fast option pricing in time-inhomogeneous Levy models us¬ 
ing partial integro differential equations (PIDEs). While large parts of the literature 
focus on numerical aspects of these pricing methods, only little is known about the 
precise link between the related deterministic equations and the corresponding con¬ 
ditional expectations representing option prices. Our main question therefore is: 
Under which conditions is there a Feynman-Kac formula linking option prices given 
by conditional expectations with solutions to evolution equations? 

In order to further specify the problem, we focus on time-inhomogeneous Levy 
models and options whose path dependency may be expressed by a killing rate. In 
this setting with A = (4()[o.t] the Kolmogorov operator of a time-inhomogeneous 
Levy process, killing rate (or potential) k : [0, TjxR* 1 ^ R, source f : [0, T]xR d -> 
R and initial condition g : R d —> R, the Kolmogorov equation is of the form 

d t u + A r-tu + n T -tU = f, 

u(0)=g. ^ 

Adopting a heuristic approach, one would typically assume that equation ([TJ) has 
a classical solution u. If this solution is sufficiently regular to allow for an application 
of Ito’s formula and moreover satisfies an appropriate integrability condition, the 
following Feynman-Kac-type representation 

T 

u{T-t,L t ) = E^g{L T )e-^ KhiLh)dh + Jf(T-s,L s )e~ f‘ KhiLh)dh ds^F^ (2) 

t 

follows by standard arguments and taking conditional expectations, see equations 
m and ED on page l30l for a detailed derivation. Then, the conditional expectation 
© can be obtained by solving Kolmogorov equation © by means of a deterministic 
numerical scheme. Such an argumentation hinges on a strong regularity assumption 
on the solution u and thus implicitly on the data of the equation, g, /, A and 


FEYNMAN-KAC FORMULA FOR LEVY PROCESSES 


3 


k. We have to realize, however, that this constitutes a serious restriction on the 
applicability of such a heuristic approach. 

To do justice to the complexities of financial applications, we pay special atten¬ 
tion to identifying appropriate conditions for the validity of equation m for financial 
applications. Often, discontinuous killing rates constitute a natural choice, as we 
will show in several detailed examples in Section [5] In particular indicator functions 
as killing rates turn out to be key to a wide variety of applications, both in math¬ 
ematical finance and in probability theory. As one typical application we propose 
and study a flexible family of employee options in Section 15.11 and illustrate the 
numerical effect of such killing rates in Section [6] The fundamental role of killing 
rates of indicator type is killing the process outside a specified domain, which makes 
them attractive for applications. Moreover, they are closely related to occupation 
times and exit times of stochastic processes as we outline in Sections 15.31 and 15.41 
We furthermore find that discontinuous killing rates form a common root of exit 
probabilities of stochastic processes and the distribution of supremum processes. 
As such they apply to the prices of path-dependent options like those of barrier, 
lookback, and American type. In view of these considerations, which are both of a 
theoretical and applied nature, we will also want to allow for non-smooth and even 
discontinuous killing rates in Kolmogorov equation (111- 

Discontinuities in the killing rate n result in non-smoothness of the solution u 
of Kolmogorov equation 0- In particular, one cannot expect u £ C 1,2 . Assume 
u( 0) ^ 0 and n = l(_ O o,0)‘ i ©i then x >->■ u(t,x) £ C 2 implies x n(x)u(t, x) £ 
C, which obviously is a contradiction. Hence, for our purposes, the assumption 
that Ito’s formula can be applied to the solution u is futile. Neither is it reasonable 
to assume that equation © has a classical solution. Let us emphasize that such 
irregularity is not only inherent in equation © if the killing rate is discontinuous, 
but also a typical feature of Kolmogorov equations for other path-dependent option 
prices. Prominent examples are boundary value problems related to barrier options 
in Levy models as well as free boundary value problems for American option prices. 
In each of these cases, the use of a generalized solution concept is called for. 

Among the possible generalizations of classical solutions of partial differential 
equations, we find that viscosity and weak solutions are the ones that are most 
commonly discussed. Viscosity solutions directly abstract from pointwise solutions 
by introducing comparison functions that are sufficiently regular, while the root of 
weak solutions is the problem formulation in a Hilbert space. Conceptually, both 
have their advantages. From a numerical perspective, viscosity solutions relate to 
finite difference schemes, whereas weak solutions are the theoretical foundation of 
Galerkin methods, a rich class of versatile numerical methods to solve partial dif¬ 
ferential equations. Relying on their elegant Hilbert space formulation, Galerkin 
methods by their very construction lead to convergent schemes as well as to a lucid 
error analysis. They furthermore distinguish themselves by their enormous flexi¬ 
bility towards problem types as well as compression techniques. Both theory and 
implementation of Galerkin methods have experienced a tremendous advancement 
over the past fifty years. They have become indispensable for today’s technological 
developments in such diverse areas as aeronautical, biomechanical, and automotive 
engineering. 

In mathematical finance, Galerkin pricing algorithms have been developed for 
various applications, even for basket options in jump models. Furthermore, numer¬ 
ical experiments and error estimates have confirmed their efficiency both in theory 
as well as in practice. See (2(1), and e.g. Matache, von Petersdorff and Schwab 
(2004), Matache, Schwab and Wilder (2005), von Petersdorff and Schwab (2004). 
We present the implementation of a related Galerkin method to price call options 


4 


K. GLAU 


adjusted with a killing rate in Section [G] Furthermore, Galerkin-based model re¬ 
duction techniques have a great potential in financal applications, see Cont, Lantos 
and Pironneau (2011), (33), and (41), Haasdonk, Salomon and Wohlmuth (2012) 
and (l23l l. 

Feynman-Kac representations for viscosity solutions with application to option 
pricing in Levy models have been derived in ©) and (fl4h . Results linking jump 
processes with Brownian part to variational solutions had already been proven ear¬ 
lier in (Q). However, in order to cover some of the most relevant financial models, 
we have to consider pure jump processes, i.e. processes without a Brownian com¬ 
ponent, as well. Pure jump Levy models have been shown to fit market data with 
high accuracy and have enjoyed considerable popularity, see for instance a, 
(12). Moreover, statistical analysis of high-frequency data supports the choice of 
pure jump models, see (jT|) . 

We realize that pure jump processes differ significantly from processes with a 
Brownian part. The Brownian component translates to a second order derivative 
in the Kolmogorov operator, while the pure jump part corresponds to an integro 
differential operator of a lower order of differentiation. Accordingly, the second or¬ 
der derivative is only present in Kolmogorov operators of processes with a Brownian 
component. As a consequence, the solution to the Kolmogorov equation of a pure 
jump Levy process does not lie in the Sobolev space if 1 , the space of quadratic 
integrable functions with a square integrable weak derivative. Therefore we need 
a more general solution space. In order to make an appropriate choice, recall that 
Levy processes are nicely characterized through the Levy-Khinchine formula by the 
Fourier transform of their distribution or, equivalently, by the symbol. Moreover, 
the symbol is typically available in terms of an explicit parametric function and as 
such is the key quantity to parametric Levy models. For a wide range of processes, 
the asymptotic behaviour of the symbol ensures that the solution of the Kolmogorov 
equation belongs to a Sobolev-Slobodeckii space, i.e. it has a derivative of fractional 
order. Even more, parabolicity with respect to Sobolev-Slobodeckii spaces of the 
Kolmogorov equations related to Levy processes has been characterized in terms of 
growth conditions on the symbol in (121). 

So as to allow for typical initial conditions, such as the payout function of a 
call option in logarithmic variables and the Heaviside step function that relates 
to distribution functions, we base our analysis more generally on exponentially 
weighted Sobolev-Slobodeckii spaces. We therefore generalize the characteriza¬ 
tion of parabolicity to time-inhomogeneous Levy processes and to exponentially 
weighted Sobolev-Slobodeckii spaces. In HI) existence and uniqueness of weak so¬ 
lutions in exponentially weighted Sobolev-Slobodeckii spaces of Kolmogorov equa¬ 
tions related to time-inhomogeneous Levy processes and a Feynman-Kac formula 
has been established. Here, we generalize these results to solutions of Kolmogorov 
equations related to time-inhomogeneous Levy processes with possibly discontinu¬ 
ous killing rates. Technically, the present setting is more difficult since the Fourier 
transform of the solution is not explicitly available and, moreover, the solution is 
not sufficiently regular for an application of Ito’s formula. 

The fruitful relation between pseudo differential operators (PDOs) and Markov 
processes via their symbols has already been extensively used to establish existence 
of stochastic processes, see for instance the monographs of Jacob from (2001), (2002) 
and (2005). For a short overview on the different approaches to construct Feller 
processes and the use of pseudo differential calculus in this context see Chapter 
III in the monograph of Bottcher, Schilling and Wang (2013). Let us observe that 
our question is of a different nature: We establish a Feynman-Kac-type represen¬ 
tation of the form m, while existence of the stochastic processes involved, L and 
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the conditional expectation, are known. An interesting feature of our approach is 
that we do need not impose growth conditions on the (higher-order) derivatives 
of the symbol as in the standard symbolic calculus. Our approach is more closely 
related to (1271 ) , where a class of martingale problems is solved tracing back the ex¬ 
istence of the processes to parabolicity of the Kolmogorov equations with respect to 
anisitropic Sobolev-Slobodeckii spaces. Compared to the setting in ((27f ), we restrict 
ourselves to isotropic spaces and constant coefficients, but, more generally, allow 
for exponentially weighted spaces and possibly discontinuous killing rates. 

To comprise all of the requirements, we state our research question more precisely 
as follows: Under which conditions on the time-inhomogeneous Levy process L, the 
possibly discontinuous killing rate k, the source f and initial condition g is there 
a unique weak solution in an exponentially weighted Sobolev-Slobodeckii space of 
Kolmogorov equation © that allows for a stochastic representation of form m? 

To answer our research question, we introduce in the next section the necessary 
notation and concepts. We use this framework first to characterize parabolicity 
of the Kolmogorov equation in terms of properties of the symbol in Theorem 13.31 
Prepared thus, we formulate our main result, the Feynman-Kac-type representation 
of the weak solution of Kolmogorov equation © in Theorem 13.41 In Section 0] we 
find that it is a wide and interesting class of stochastic processes that fall within the 
scope of this result. Analysing its applications in Section [5] leads us from typical 
financial problems further to the characterization of purely probabilistic objects and 
finally back to the original quantum mechanical ideas of Feynman and Kac—yet in 
a relativistic guise. Exploiting the advantages of Theorem 13.41 further, we return 
to its practical realization and implement a Galerkin scheme to solve Kolmogorov 
equation © in Section[G] We find that thanks to Theorem l3.4l the solutions obtained 
thus correspond to option prices. With the numerical implementation at hand, we 
visualize and discuss the effect of killing rates of indicator type. Section [7] presents 
a robustness result for weak solutions that is required in our proof of Theorem 13.41 
in Section [8] In this last section we also identify desirable regularity properties 
of the solutions to the Kolmogorov equation. Appendix [X] provides two technical 
lemmata for the symbol and the operator, and Appendix [B] concludes with the 
proof of Theorem 13.31 


2. Preliminaries and notation 

In order to present the main result of the present article, we first introduce the 
underlying stochastic processes, the Kolmogorov equation with killing rate, its weak 
formulation as well as the solution spaces of our choice. We denote by Cg° (R d ) the 
set of smooth real-valued functions with compact support in R d and let 

T(ip) := ( ip(x) dx (3) 

R d 

be the Fourier transform of <p G Co°(R d ) and J 7-1 be its inverse. 

Let a stochastic basis (H, J 7 , {J-’t)o<t<T,P) be given and let L be an Revalued 
time-inhomogeneous Levy process with characteristics (b t ,at,F t ;h)t> o- That is L 
has independent increments and for fixed t > 0 its characteristic function is given 

by 

= e~ti As( " i « )ds for every f G R d , (4) 

where, for every t > 0 and £ G R d , the symbol of the process is defined as 

MO ■= ^(OM) +i(Z,bt) - J (e~ l{i ’ v) -1 + *(£,%)}) F t (dy). 

R d 


( 5 ) 
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Here, for every s > 0, a s is a symmetric, positive semi-definite dxd-matrix, b s £ R d , 
and F s is a Levy measure, i.e. a positive Borel measure on R d with F s ({0}) = 0 
and J R d(\x \ 2 A l)F s (dx) < oo. Moreover, h is a truncation function i.e. h : R d — s- R. 
such that h(x)F t (dx) < oo with h(x) = x in a neighbourhood of 0. We 

assume the maps s H> er s , s >->• b s and si-> J {\x\ 2 A l)F s (dx) to be Borel-measurable 
with, for every T > 0, 

T 

J (|&«| + IWI-M(dxd) + J(\x\ 2 Al)F s (dx)j ds < oo, (6) 

0 Rd 

where || • || M(dxd) is a norm on the vector space formed by the d x d-matrices. 

The Kolmogorov operator of the process L is given by 


A t <p(x) 



(7) 


for every ip £ Cq°( R d ), where hj denotes the j-th component of the truncation 
function h. By some elementary manipulations we obtain 

Atip = F~ l (A t F(g>)) for all tp £ C 0 °°(R d ), (8) 


which shows us that the Kolmogorov operator A is a pseudo differential operator 
with symbol A. 

Following the classical way to define solution spaces of parabolic evolution equa¬ 
tions, we introduce a Gelfand triplet (V, F [, V*), which consists of a pair of separable 
Hilbert spaces V and F[ and the dual space V* of V such that there exists a con¬ 
tinuous embedding from V into H. We then denote by L 2 (0,T;.ff) the space of 

weakly measurable functions u : [0, T] —> FI with ||u(i)||jydi < oo and by d t u 
the derivative of u with respect to time in the distributional sense. The Sobolev 
space 


W\0,T-,V,H) :={u£ L 2 (0 ,T; V) | d t u £ L 2 (0 ,T; W)}, 


(9) 


will serve as solution space for equation ©■ For a more detailed introduction to 
the space W 1 (0, T; V, iL), which relies on the Bochner integral, we refer to Section 
24.2 in ([47). More information on Gelfand triplets can be found for instance in 
Section 17.1 in (l47h . 

Usually, variational equations of a similar type as the heat equation are for¬ 
mulated with respect to Sobolev spaces, and thus are based on both H 1 and L 2 . 
Since we include pure jump processes in our analysis, operator 0 may be of frac¬ 
tional order. We therefore work with Sobolev-Slobodeckii spaces, which formalize 
the notion of a derivative of fractional order. Turning to a typical financial prob¬ 
lem, we express the price of a call option as solution to a Kolmogorov equation of 
type ©■ We then obtain k = 0 and / = 0, while the initial condition is given by 
g(x) = (Sq e x —K) + . We now have to realize that the initial condition g £ L 2 (R d ) 
and we cannot use an L-based approach. The exponentially dampened function, 
x 1 —^ g(x) e vx , though belongs to L 2 ( R d ) for every g < —1. Thus, in order to incor¬ 
porate initial conditions that typically arise in financial problems, we allow for an 
exponential weight. We further increase the class of function spaces by a domain 
splitting argument, see Remark [5© on pagelTTl 
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To make these considerations formally precise, we define the exponentially weighted 
Sobolev-Slobodeckii space Hf(R d ) with index a > 0 and weight 77 £ R d as the com¬ 
pletion of Co°(R d ) with respect to the norm || ■ || h™ given by 

m 1 % ■■= / (i+ieirwtt-^rde. (10) 

R d 

Observe that this is a separable Hilbert space. For rj = 0 the space HTjR d ) coincides 
with the Sobolev-Slobodeckii space H a (R d ) as it is defined e.g. in (1471 1. For a = 0 
the space H“( R d ) coincides with the weighted space of square integrable functions 
L 2 (R d ) := { u £ Tjoc^lR/*) | x i-A u(x) e l ' 71 ’ x ' 1 £ L 2 (R d )}. Furthermore, we denote the 
dual space of iL“(R d ) by 

Let a : [0, T] xff“(R d ) xiL“(R d ) —> R be a family (a t ) t£ [o l r] of bilinear forms that 
are measurable in t with associated linear operators At ■ H“(R d ) — > (lL“(R d ))* 
given by 

At{u)(v) = a t {u, v) for all u,v £ ff“(R d ) ( 11 ) 

and whose related symbols A t : R d —> € are such that 

M<fi) = J A t (C)^)(C)dC for all ^ £ Q°(R d ). ( 12 ) 

R d 

We close the section with the weak formulation of Kolmogorov equation 0 - 

Definition 2.1. LetV = R d ) and H = L 2 (R d ), k : [0,T]xR d —> R measurable 
and bounded, f £ L 2 (0,T;V*) and g £ H. Then u £ W 1 (0, T; V, H) is a weak 
solution of Kolmogorov equation 0, if for almost every t £ (0, T), 

(d t u(t),v) H + a T -t(u(t),v) + {nT-tu(f),v) H = (f(t) \v)vxv forallv£V (13) 
and u(t ) converges to g for t f 0 in the norm of H. 

3. Main results 

Equipped with the necessary notation and concepts, we now focus on our main 
purpose, providing a Feynman-Kac formula linking weak solutions of PIDEs with 
killing rates to conditional expectations. 

Following a classical way to prove existence and uniqueness of weak solutions 
of a parabolic equation, we verify continuity and a Garding inequality of its bilin¬ 
ear form. We specify the notion of parabolicity accordingly and adapt it to our 
framework: 

Definition 3.1. Let A be an operator associated with bilinear form a. 

We say A, respectively a, is parabolic with respect to H^ 2 (R d ),L 2 (R d ), if 
there exist constants C, G > 0, G' > 0 such that uniformly for all t £ [0, T] and all 
u,v £ H% /2 { R d ), 

|a t (w,u)| < C\\u\\ H c./ 2 ( Rd) \\v\\ H «/ 2 ( Rd) (Continuity (Cont-a)) 

a t (u, u) > G|M|^/ 2(Rd) - G'\\u\\ 2 L ^ nd) . (Garding inequality (Gard-a)) 

For R C R d , we say that the parabolicity of A, respectively a, (with respect to 
(H^ 2 (R d ), L 2 (R d )) veR ) is uniform in [0,T] x R, if for all u,v £ U V £rH^ ~ (R d ) 
the mapping t at{u,v) is cadlag and there exist constants C,G > 0, G' > 0 
such that uniformly for all rj £ R, all t £ [0,T] and u,v £ H^ 2 (R d ) inequalities 
(Cont-a) and (Gard-a) are satisfied. 
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As highlighted in equation Q. the Kolmogorov operator of a time-inhomo¬ 
geneous Levy process is a pseudo differential operator. Its symbol is explicitly 
known for various classes and in general is characterized by the exponent of the 
Levy-Khinchine representation. Therefore, we express our main assumptions in 
terms of the symbol of the process. For Levy processes with symbols A , it has been 
shown in (2ll). Theorem 3.1, that the corresponding bilinear form is parabolic with 
respect to H a / 2 (R d ), L 2 (R d ) if and only if constants C,G,G' > 0 and 0 < /3 < a 
exist such that for every £ £ lR, d , 

|A(£)| <C7(l + |£|)“ ( 14 ) 

>G(l + |C|) a -G'(l + |^|) /3 . (15) 

We generalize this growth condition so as to render it suitable for the setting of 
time-inhomogeneous Levy processes and weighted Sobolev-Slobodeckii spaces. We 
find that an extension of the bilinear form to weighted Sobolev-Slobodeckii spaces 
corresponds to a shift of the symbol in the complex plane. Symbols can be extended 
to complex domains if the appropriate exponential moment condition is satisfied. 
Let L be a time-inhomogeneous Levy process. First notice that L t is infinitely 
divisible with Levy measure F t (dx) := f* F s (dx) ds for every t £ [0, T], as has been 
shown by Lemma 1. Theorem 25.17 in (f43h now implies that, for all 77 E R d , 


T 


// 

0 |x|>l 


eh)*®) F)(da;) dt < oo 


{EM {if}) 


is equivalent to the exponential moment condition ] < oo and 

E [ e (H+ri,L t ) ] = e -/ 0 *^(-C+iu)d« for aU £ e & nd t > 0. (16) 

We therefore formulate the conditions in terms of an exponential moment condition 
on the process and growth conditions on the symbol extended to a complex domain. 
It turns out that this complex domain can conveniently be chosen as a tensorized 
complex strip. More precisely, for weight 77 = ( 771 ,..., rid), let 

U v := {z £ <D d | S(zj) £ {0} U sgn(? 7 j)[ 0 , |) for j = (17) 

R v := sgn(? 7 i)[ 0 , |r?i|] x ■ • ■ x sgn(? 7 d )[ 0 , \r] d \\. (18) 

From Theorem 25.17 in (43), we also know that the complex set on which A s is 
definable is convex. Lemma 2.1 (c) in Hi shows for the present setting that the 
map z 1 —y A t (z) has a continuous extension to the complex domain U- v that is 

O 

analytic in the interior 17-^. 

We will derive the main results related to the following set of conditions. 

Conditions 3.2. For weight 77 £ !R d and index a £ (0,2], let A = (A t ) tg [ 0 ' 7 ’] be a 
symbol with extension to U- v and, if available, let L denote the time-inhomogeneous 
Levy process with symbol A. 

(Al) For every 77 ' £ R rn 

F[e - ^ ,Lt ' > ] < 00 . (Exponential moment condition {EM)) 

(A2) There exists a constant C > 0 such that uniformly for all 77 ' £ R^ and 
t e [0,T], 


A t ^- W )\ <^(1 + 1 ^ 1 )“. 


(Continuity condition (Cont-A)) 
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(A3) There exist constants G > 0, G' > 0 and 0 < /3 < a such that uniformly for 
all rf G R v and t G [0, T], 

!R(At(£ - it]')) > G(l + |(|) Q - G'(l + Ifl)^. (Garding condition (Gard-A)) 

(A4) For every fixed r/ € R v and £ G R d the mapping 1 1 -> A t (£ — irf ) is cadlag. 

We say that A has Sobolev index a uniformly in [0, T] x R v , if A has an extension 
to U- v that satisfies (A2)-(A4). If A is the symbol of process L , we also say L has 
Sobolev index a uniformly in [0,T] x R v . 

Conditions (Al)-(A4) are satisfied for a large set of processes, for instance 
for tempered stable and normal inverse Gaussian processes as well as their time- 
inhomogeneous extensions. In Section [4] we look in detail at the verification of the 
conditions. 

Notice that for rj = 0 we have R v = {0}. Thus (Al) is trivially satisfied and (A2)- 
(A3) simplify accordingly. This case corresponds to the case of Sobolev-Slobodeckii 
spaces without weighting and is covered by the following results. If, moreover, the 
symbol is constant in time, (A4) is trivially satisfied and (A1)-(A4) reduce to (fT4l) 
and (USD- Conditions (Al)-(A3) were introduced in (|l(fl to show existence and 
uniqueness of weak solutions of the related Kolmogorov equation (without killing 
rate) along with a Feynman-Kac-type formula with application to European option 
prices in time-inhomogeneous Levy models. We additionally require (A4), which 
only imposes a mild technical restriction. 

Our framework defined, let us now state our main results. We first show in 
Theorem l3.3l the equivalence between parabolicity with respect to weighted Sobolev- 
Slobodeckii spaces and growth conditions (A2) and (A3), thereby generalizing the 
result for Sobolev-Slobodeckii spaces and conditions and m to the present 
setting. The characterization of (uniform) parabolicity in terms of conditions on 
the symbol is interesting in its own right. It is, moreover, one of the key steps in our 
proof of Theorem l3.4l below. which establishes a Feynman-Kac-type representation. 

Theorem 3.3. For r) G R rf and a G (0,2], let L be a time-inhomogeneous Levy 
process satisfying exponential moment condition (Al). Then the following two as¬ 
sertions are equivalent. 

(i) The Kolmogorov operator of L is uniformly parabolic in [0, T] x R v with respect 

(ii) L has Sobolev index 2 a uniformly in [0,T] x R v . 

The proof of Theorem l3.3l is given in AppendixjBj where, moreover, Theorem lB.il 
provides a more general version of this result for operators and symbols that are 
not necessarily related to stochastic processes. 

Assume (A1)-(A4). Then Theorem 13.31 shows in particular the parabolicity of 
the related bilinear form uniformly in time with respect to H(f( R d ) and L^(R d ). 
Now the classical existence and uniqueness result, see for instance Theorem 23.A 
in ©, gives us that Kolmogorov equation © has a unique weak solution u in the 
space W 1 (0,T;iL“(R d ), L*(R d )). 

We now turn to the stochastic representation of this solution. For an integrable 
or nonnegative random variable A' we denote 

Eo, x (X) := E x (X), E tjX (X) := E(X\L t = x) for t > 0, (19) 

where x i—>■ E(X\L t = x) is the factorization of the conditional expectation E(X\L t ) 
and E x the expectation with respect to the probability measure P x such that 
P X {L 0 = x) = l. 
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Theorem 3.4. For g £ R d and a £ (0,2], let L be an M. d -valued time-inhomo- 
geneous Levy process with symbol A = (^4t)tg[o,T] that satisfies (Al)-(Af). Then 

( i ) for k : [0,T] x WL d —> R measurable and bounded, f £ L 2 (0,T-(H^ 2 (R d ))*) 
and g £ L 2 (H d ) Kolmogorov equation |1]) has a unique weak solution u £ 
W 1 (0, T; Hy/ 2 (H d ), L 2 (R d )'j; 

(ii) if, additionally, f £ L 2 (0, T; for some l > 0 with l > (d — a)/2, then, 

for every t £ [0,T] and a.e. x £ R d , 

u(T -t,x)= E t , x (g(L T ) e" & Kh{Lh) dh 

T r \ ( 20 ) 

+ / f(T — s, L s ) e~ ht K h{ L h) dh d s j _ 


As we have already seen, part (i) of Theorem 13.41 follows from Theorem 13.31 and 
the classical existence and uniqueness result for solutions of parabolic equations. 
Part (ii) is considerably more involved and Section [5] is devoted to its proof. 

The major benefit of Theorem 13.41 for financial applications is that conditional 
expectations of form GOD, which naturally appear as derivatives and asset prices, 
are now characterized by weak solutions of PIDEs. Therefore, the prices can be 
computed by numerically solving an equation of form ©■ So as to illustrate the 
method and the effect of killing rates, we present among others an application to 
employee options in Section [5] and provide its Galerkin discretization in Section [G] 
Theorem l3. 41 furthermore shows a specific type of regularity of conditional expec¬ 
tation GOD- It is interesting to identify sufficient conditions for Holder continuity. 
Theorem 8.2 in Nezza, Palatucci and Valdinoci (2011) provides the appropriate 
Sobolev embedding result. Thus we obtain as an immediate consequence of Theo¬ 
rem [331 the following corollary. 


Corollary 3.5. Under the assumptions and notations of Theorem in the uni¬ 
variate case, i.e. for d = 1, for a £ (1,2] and any fixed t £ (0,T), the function 
x i ^ u(t,x) is X-Holder continuous with A = ^-=2. , i.e. 


\u(t, x) — u(t, y)\ 
sup . 

i,v€R,iys \ x ~ y I 


< 00 . 


In particular, x u(t,x) is continuous and equality (BUI) in Theorem \3.4\ holds for 
every x £ R. 


4. Examples of classes of time-inhomogeneous Levy processes 

Let us explore the nature of Conditions (Al)-(A4) and show that they are sat¬ 
isfied for a wide class of processes. Conditions (Al)-(A4) naturally apply to pro¬ 
cesses that are specified through their symbol. Notice that the symbol is expressed 
in terms of the characteristics of the process. We exploit this in Proposition 14.71 
to establish concrete accessible conditions for real valued time-inhomogeneous pure 
jump Levy processes with absolutely continuous Levy measures, while Proposition 
14.31 treats time-inhomogeneous multivariate jump diffusions. 

We should, however, realize that Conditions (Al)-(A4) are not satisfied by all 
Levy processes. On the one hand, continuity and Garding condition (A2) and (A3) 
have implications for the distributional properties of the process: 

Remark 4.1. Fix rj £ R d and a £ (0, 2], and let L be a time-inhomogeneous Levy 
process with symbol A = (A t )t> o- If Garding condition (AS) is satisfied for weight 
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p and index a, then there exist C\, Gi > 0 such that uniformly for all rf £ R v and 
0 <s<t<T, 


g- Is A u (£-i7]')du I < ^ e -(t-s)C 2 |?r _ 


( 21 ) 


In particular, (AS) implies for every t € (0,T] that the distribution of L t has a 
smooth Lebesgue density. 

On the other hand, continuity and Garding condition (A2) and (A3) relate to 
the path behaviour of the process: If a Levy process with symbol A satisfies (A2) 
and (A3) for a £ (0, 2) and p = 0, then a is its Blumenthal-Getoor index, as shown 
in (|2L). Theorem 4.1. Hence, every pure jump Levy process satisfying assumptions 
(A2) and (A3) has infinite jump activity. On this basis we may for instance conclude 
that compound Poisson processes do not satisfy (A3). 

Variance Gamma processes have Blumenthal-Getoor index 0 and thus do not 
satisfy both (A2) and (A3), as noticed in part (iv) of Example 4.1 in ©)■ However, 
pure jump Levy processes can be approximated by a sequence of Levy jump diffusion 
processes with nonzero Brownian part. This can always be achieved by adding a 
diffusion part and letting its volatility coefficient tend to zero. Example 14.41 shows 
that pure jump Levy processes can be approximated by Levy processes for which 
(A1)-(A3) are satisfied for weight p = 0 and index a = 2. (JJ) provide a sequence 
with better approximation properties for Monte Carlo techniques, which could be 
exploited further. 

Before continuing the discussion on the validity of Conditions (Al)-(A4) for 
other classes of processes, we observe the following. 

Remark 4.2. For p £ lR, d and a £ (0, 2] , let A be the symbol of a Levy process that 
satisfies exponential moment condition (Al). By virtue of Lemma \ A. 1\ in Appendix 
HI and the continuity of symbols of Levy processes (as mappings from R d to C), the 
validity of continuity condition (A2) for A is equivalent to the following asymptotic 
condition: For every N > 0 there exist a constant G > 0 such that for every p' £ R v , 


3i(A(£ — ip')') > G|£|“ — A(ip') for every £ £ R d such that |£| > N. (22) 


We devote the remainder of this section to providing sufficient conditions for the 
validity of Conditions (Al)-(A4) for time-inhomogeneous jump diffusions, for pure 
jump Levy processes, and for time-inhomogeneous processes. 


4.1. Jump diffusions. For time-inhomogeneous Levy jump diffusion processes we 
find that Conditions (A1)-(A4) are satisfied under remarkably weak conditions: 

Proposition 4.3. Fix some p £ III d . Let L be a time-inhomogeneous Levy process 
with characteristics (b t , at, F t ', h)o<t<T such that 


sup 

te[o,T] 



for every p' £ R v and (23) 


M>i 



(24) 


Then (A1)-(A3) are satisfied for weight p and index a = 2. 
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Proof. Due to the equivalence of EM(y) and the exponential moment condition, 
(121 implies (Al). Observing that 

- in')) ={bt,v') + + J (((^OW) - l)e" <I ''’ x> -l)F t (dx) 

n d 

+ ^<£>0*0 + J (cos((£, ir)) —ljF t (dx) 

M d 

and f Rd (cos((£, a;)) — l)F t (dx) > 0, inequalities (121 and (121 yield Garding con¬ 
dition (A3) with a = 2. Similarly, inequalities (E3l) . (HH1) and 

Q(A t ({ - ir/')) =(bf v ,£)+ f (sin ((£,x)) - (£,h(x))) e~^' x) F t (da;), 

where bf v = b + a ■ y' + f Rd (e^ ,v ^ — l)h(y)F(dy) as defined in Lemma lA.ll in 
Appendix iBl yield continuity condition (A2), which concludes the proof. □ □ 

For Levy jump diffusion processes the conditions simplify considerably: 

Example 4.4 (Multivariate Levy processes with Brownian part). Fix y £ R d and 
let L be an R d -valued Levy processes with characteristics {b,a,F;h) such that a is 
a positive definite matrix and the Levy measure F satisfies Jj a .| >1 e~ vx F t (dx) < oo. 
Then (Al)-(A3) hold for weight y £ R d and index a = 2. 

In order to verify those assumptions of Proposition 14.31 that concern the pure 
jump part of the process, it suffices to consider the pure jump processes separately, 
as the following Lemma shows. 

Lemma 4.5. For j = 1,2, let L 3 be two stochastically independent time-inhomo- 
geneous Levy processes with symbol A 3 such that (Al)-(Af) are satisfied for the 
same weight y £ R d and the possibly different indices a 3 . Then the sum L := L 1 +L 2 
is a time-inhomogeneous Levy process with symbol A := A 1 + A 2 , and (Al)-(Af) 
are satisfied for weight y and index a := max(ar,ar). 

Lemma [H generalizes Remark 4.1. in (0) to the case where y 0 and we omit 
its elementary proof. 

4.2. Pure jump Levy processes and operators of fractional order. We now 

consider a class of multivariate processes, which frequently occurs in finance and 
whose symbol is explicitly given. 

Example 4.6 (Multivariate Normal Inverse Gaussian (NIG) processes). Let L 
be an 1 R, d -valued NIG-process, i.e. a Levy process such that L\ = (L\,...,L d ) ~ 
NIGd(u, /3, S , //, A) for parameters a, S > 0, /3, p £ R d and symmetric positive defi¬ 
nite matrix A £ R dxd with a 2 > (/3, A/3). The symbol of L is given by 

A{u) = i(u, n) - s(^\/a 2 — ((3, A/3) — \J a 2 - {(3 + iu, A (/3 + iu))j , 

where we denote by (■, ■) the product (z,z r ) = )Cj=i z j z j f or z £ Compare e.g. 
equation (2.3) in \2A) . 

Assumptions (A1)-(A3) are satisfied for index a = 1 and every y £ R d such that 
a 2 > (/3 + y', A(/3 + y’)) for all y' £ R v . This is in particular the case, if 

||^|| 2 + IH| 2 <d 2 /l|A||. (25) 

To summarize, if the parameters of L satisfy m, Conditions (Al)-(Af) are sat¬ 
isfied for weight y and Sobolev index 1 . 





FEYNMAN-KAC FORMULA FOR LEVY PROCESSES 


13 


Since pure jump Levy processes can be defined through a Levy measure and 
a constant drift, we are interested in finding conditions on both the Levy mea¬ 
sure and the drift that imply Conditions (A1)-(A4). Let us address this issue for 
real-valued time-homogeneous pure jump Levy processes whose Levy measure is 
absolutely continuous. For this class we generalize Proposition 4.2 in (21) to time- 
inhomogenuity and weights i] ^ 0. Thereby we obtain explicit conditions on the 
characteristics that imply Conditions (Al)-(A4). 


Conditions 4.7. Fix rj £ ]R and a £ (0,2], and let L be a real-valued time- 
inhomogeneous Levy process with characteristics (bt, at, F t ] h) t > 0 . 

( F1 ) Jo J|x|>i e ~ VX F ‘(dx) dt < °°» 

(F2) F t is absolutely continuous for every t £ [0,T] with density ft, i.e. F t (dx) = 
ft(x ) dx. Denote the symmetric part by ff ym (x) := (ft(x) + /(—x))/2 and 
the antisymmetric part by f? sym (x) := ft(x) — ff ym (x). 

(F3) There exist constants C\,C 2 ,e > 0 and 0 < /3 < a < 2 and a function 
g : [0, T] x [—e, e] —>• ]R such that uniformly for all t £ [0, T\, 

ft Vm (x) < +g(t,x) and \g(t,x)\ < for all |x| < e. 

(F4) If a = 1, there exist constants C^,e > 0 and (3 £ (0,1] such that uniformly 
for all t £ [0, T], 

\ft sym (x)\ < |-^ for all |x| < e. (26) 

If a < 1, then inequality © holds for some (3 £ [0, a] and, moreover, 
bt = J h(x)F t (dx) for every t £ [0,T]. 

Proposition 4.8. Let L be a real-valued time-inhomogeneous pure jump Levy pro¬ 
cesses with characteristics (bt , 0, F t ) t > 0 . Then, 

(i ) Condition (FI) is equivalent to (Al); 

(ii) Conditions (F1)-(F3) imply (Al) and (A3); 

(Hi) Conditions (F1)-(F4) imply (Al)-(A3). 


Proof, (i) Since d = 1, we have R v = sgn(?7)[0, r)\, and part (z) directly follows 
from Theorem 25.17 in ( |43h . 

(ii) We denote f t , v (x) := f t (x), ff ym (x) := (f t , v (x) + ft, v (~x))/ 2 and 
yasym jsym _ fsy™ Then the elementary identity ab + cd = (a + c)(b + d) /2 + 
(a — c)(b — ’d)/ 2 yields 

ft V v m ( x ) = C0S Hvx)ft V rT( x ) + sin Hvx)ft" v ym (x). 

We notice that there exists a constant c > 0 such that 


cosh(7yx) > e ve and | sinh(? 7 x)| < c|x| for every |x| < e. 


Moreover, since / > 0, the triangle inequality implies that \ff^ vm (x)\ < ff v v m (x ) 
for every (t,x) £ [0,T] x 1R.. This shows that Condition (F3) also remains valid 
when we replace ff ym by f f ym . Now part (ii) follows from inequality (4.16) in the 


proof of Proposition 4.2 in (21). 

(iii) Along the same lines as in the proof of part (ii), we observe that 


ft,7 m ( x ) = sinh(r;x)/ t s ^ m (x) + cosh(r?x)/ t “* ym (x). 

Thus, the validity of Condition (F4) also remains valid when replacing ff svm by 
fasym qq leil) Proposition 4.2 in ([2j]) shows the assertion for a = 1. 
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For a < 1 we have bt 
notation therein, bf' 1 = 


= f h(x)F t (dx). According to Lemma [AJj and using the 
/ h{x)Ffi v (da;). Hence, for all t £ [0,T] and r( £ R v , 


7m(A(£ — it]')) = j sin(£a;)e vx F t (dx). 


Estimating the real part of the A(£ — irf) along the same lines as in the proof of 
Proposition 4.2 in 421), we obtain continuity condition (A2). □ □ 


We now apply Proposition l4.8l to a concrete class of processes, which is frequently 
used to model asset prices: 


Example 4.9 (Univariate generalized tempered stable Levy process). Let L be a 
generalized tempered stable Levy process with parameters C_, C+ > 0 such that 
C- + C+ > 0 and G, M > 0 and Y_, Y + < 2. That is, A is a pure jump Levy process 
with characteristic triplet (b, 0, F temp ; h) with E temp (da;) = f temp {x)dx, where 

( C - e Gx 

/ temP (*) = 


1*1 


1*1 


T+yZ 

’C+ 

TFTI 


Mx 


for x < 0 
for x > 0. 


For C± = 0 we set Y± := 0. For C = C- = C+ and Y = Y_ = Y + this class is 
known as CGMY, after Carr, Geman, Madan and Yor. Tempered stable processes 
are also referred to as Koponen and KoBoL in the literature, see e.g. (@). For the 
general setting see for instance 0 ). 

By Proposition \4-8[ Conditions (A1)-(A3) are satisfied for weight rj £ (— G,M) 
and Sobolev index a := max{Y + , Y_} in each of the following cases: 

(■ i ) a = max{Y + , Y_} > 1, 

(m) Y := Y_ = Y+ = 1 and C_ = C+, 

{in) 0 < a = max{Y + , Y_} < 1 and b = f h{x)F{dx). 


Further examples are considered in 0) for the case rj = 0. There, Examples 
4.5-4.7 give conditions on the parameters that imply Conditions (A2) and (A3) for 
generalized student-t, Cauchy, generalized hyperbolic and stable processes. More¬ 
over, Section 4.3 in (21) provides a sufficient tail condition on the Levy measure for 
(A2) and (A3) to hold. 


4.3. Time-inhomogeneous processes. When modeling with Levy processes in 
finance we often need to consider the larger class of time-inhomogeneous Levy 
processes, because their flexibility in time leads to a considerably better fit to the 
time-evolution of data. We therefore propose two construction principles that lead 
to parametric families of time-inhomogeneous processes satisfying (A1)-(A4). 

First, we find it natural to define a family of time-inhomogeneous Levy processes 
by inserting time-dependent parameters into a given parametric class of Levy pro¬ 
cesses. For this class it turns out to be straightforward to show the following result. 

Lemma 4.10. LetV C and {A{p,-)) p ^-p a parametrized family of symbols. Fix 
some ri £ R d and some a £ (0,2). Let (A2) and (A3) be satisfied for A, uniformly 
for all p £ V. Then, if t p(t) measurable, then (A2) and (A3) are satisfied for 

A t{f) ■= A(p(t),£) for t £ [0, T] and £ £ U- v . 

If, moreover, (p, £) A(p, £) is continuous and 1 1 ->- p{t) is is cadlag, then (A t )t >o 
is the symbol of a time-inhomogeneous Levy process L' and also satisfies (Af). If 
additionally (Al) is satisfied for L, then it is also satisfied for L'. 

For p{t) we can for instance choose a vector of piecewise constant parameters, 
so as to incorporate different short-, mid- and long-term behaviour. 
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As another natural construction let us consider stochastic integrals of determin¬ 
istic functions with respect to Levy processes. Let L be an Revalued Levy process 
and / a deterministic L-integrable R nxd -valued function. Then 

Xt~f- L t := f f(s) d L s :=( W f k (s) d L k ) 

o V k=i i ) j<d 

defines an R"-valued semimartingale with deterministic characteristics. Denote by 
( b,c,F;h ) the characteristics of L. Applying standard arguments from the semi¬ 
martingale theory, we see that the characteristics (&j Y , q Y , F x \ h)t>o of X are given 

by 

b? = f(t)b+ f (h(f(t)x) - f(t)h(x))F( dx), 

R d 

cf = f(t)cf(t) tr , (27) 

F t x {B) = I t B (f{t)x)F(dx) for every B G B(. R d \ {0}). 

In particular, X is a time-inhomogeneous Levy process in the sense of our definition 
provided that integrability condition ([6]) is satisfied for its characteristics. Moreover, 
if A is the symbol of L, the symbol A x of X is given by 

A t (f) = A (f{t) tr £) + i{^,b{h,h,f)) for every £ G R d , (28) 

where b(h,h,f) := f Rd (h(f(t)x) — f(t)h(x))F(dx). This generalizes Example 7.6 
in (ilil) . where / : [0, oo) —>■ R+. 

Lemma 4.11. Let L be a Levy process that is also a special semimartingale and 
let A denote its symbol. Let f : [0,oo) —> R raxd be a measurable function such that 
there exist constants 0 < /*, /* with 

sup \\[.f{t)f(t) tr ]~ 1 \\ l/2 < fL 1 and sup \\f(t)f(t) tr \\ 1/2 <f*, (29) 

0 <t<T 0<t<T 

where || ■ || denotes the spectral norm. Then X := f ■ L is a time-inhomogeneous 
Levy process as well as a special semimartingale and its symbol is given by 

A x (0 = A(f(tYforalHeR n . 

Fix some p > 0, rj x G R d with \rj x \ < -p and some a > 0. If Ee p ^ Lt ^ < oo for 
some t > 0, then X satisfies ( EM(R_ rl x )). If additionally A satisfies (A2) and 
(A3) for every weight ?y G R d with \rj\ < p and index a > 0, then (A2) and (A3) 
hold for A x with the same index a and weight rj x . Moreover, if (Af) holds for A 
it is also satisfied by A x . 

Proof. From the assumptions it is immediate that / is integrable with respect to L 
and, hence, X is a semimartingale with characteristics of form m- As integrability 
condition © also follows directly, we see that X is a time-inhomogeneous Levy 
process. Since L is a special semimartingale, we have /| x , >1 |a:|F(da;) < oo, where 
F denotes the Levy measure of L, and (P2!)l) implies 

T 

J J \x\F t (dx) <Tf* J |i|F(da;) < oo. (30) 

0 |x|>l \x\>l/f. 

This shows that also X is a special semimartingale. Therefore we may choose both 
h and h as the identity so that b(h , h , /) = 0. From K51) we now obtain the equality 
Af (£) = A(f (t)* r £). The assertion as to the exponential moment condition (Al) 
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we obtain analogously to The assertions on (A2)-(A4) follow immediately 

from the continuity of Levy symbols and Lemma IA.1I □ □ 

5. Applications 

Having convinced ourselves that it is a wide and interesting class of stochastic 
processes for which Theorem 13.41 links conditional expectations with weak solutions 
of PIDEs, let us now explore the virtues of the result for applications. Starting 
with pricing problems in finance, where discontinuous killing rates arise naturally, 
we furthermore find that indicator type killing rates also help us to characterize in¬ 
teresting probabilistic objects. In all of these applications the driving process L can 
be chosen freely and we may employ jump-diffusions or pure jump processes. The 
latter are intensely used in finance. Examples are NIG and generalized tempered 
stable processes, which we have shown satisfy the assumptions of Theorem [XU Fi¬ 
nally we encounter the original ideas of Feynman and Kac in a relativistic guise. In 
the context of the relativistic Sclrrodinger equation we shall see the family of NIG 
processes in a fundamental role. 

5.1. Employee options. We propose a class of employee options that flexibly 
reward the management board according to the performance of the corporation’s 
stock price. Financial instruments used in this context are called employee stock 
options and often are based on European call options. Thus the reward depends on 
the level of the stock at specific points in time. Shareholders though typically are 
interested in the performance of the stock during the whole period. They mean to 
support management decisions that push the stock price constantly to a high level. 
Moreover, it is arguably fairer to choose the reward according to the performance 
of the stock value as relative to the market evolution. 

To make this formally precise, denote by S the d-dimensional stochastic process 
that models the stock of the company and d — 1 reference assets. Let G : —> R 

be a payout profile and k : [0, T] x R d -A- R a reward rate function. For n < 0 the 
reward turns into a penalty. Moreover, we include a continuously paid salary by 
the salary function f : [0,T] x R d —R. At maturity T the employee obtains the 
payout 

G(St) e^° T Kh( ' Sh ' ) dh , (31) 

in addition to the salary 

f(t,S t )eti KhiSh)dh dt, (32) 

which is paid at each instant t £ [0,T]. Thus, the payout profile G may depend 
on the level of the stock and the reference assets. The reward rate and the salary 
may additionally be time-dependent. Note that our analysis allows us to incorpo¬ 
rate discontinuities in the reward rate. Thus threshold and indicator type reward 
functions are allowed, which is a natural choice. Indicator type killing rates for 
instance play the role of instantaneous rewards or penalties for stock price levels in 
a specified domain. 

We further use the following notation. For x = (aq,..., Xd) £ R d , let e x := 
(e Xl ,... ,e Xd ), G{x) := G(e x ), k(-,x) := -n(-,e x ) and f(-,x ) := f(T — -,e x ). We 
assume the interest rate ( rt)t>o to be deterministic, measurable and bounded. We 
model S = (Sq e Ll ,..., Sq e L ) by a time-inhomogeneous Levy process L with local 
characteristics (5, c, F; h ) such that the no-arbitrage condition, 

b\ = r t - - J (e Xi - 1 - hi(x))F t (dx) for every i = 1,..., d, (33) 

is satisfied, where hi is the i-tli component of the truncation function h. 
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The following assertion shows that the fair price of the employee option specified 
by (EH) and (IH^ll can be computed by solving the related Kolmogorov PIDE. The 
result is an immediate consequence of Theorem 13.41 

Corollary 5.1. Let p £ R d and a £ (0,2] such that G £ L 2 (R d ) and assume 
the time-inhomogeneous Levy process L satisfies (1331) and Conditions (Al)-(Af). 
Denoting x = log(S’o), the fair price 

T 

u(T,x) := E x (&(L t ) e^ (iih{Lh) - rh) dh + J f(T - s, L s ) dh ds) 

0 

of the employee option with payout profile ©, © is given by the unique weak 
solution u £ W 1 (0 1 T-H* /2 {^ d ) 1 L 2 ri {W 1 )) of 

u + Ar-tu + Kr-tu = -/, u(0) = G. (34) 

See Section O for a numerical implementation of equation (IH-fll . 

5.2. Levy-driven short rate models. Levy driven term structure models were 
introduced first in ( 0 ). Here, we consider a short rate of the form 

r t :=r(f,L t ) (35) 

with an Revalued time-inhomogeneous Levy process L and a measurable and 
bounded interest rate function r : [0,T] x R d —> R. We allow for discontinuities in 
the function r and thus for thresholds in factor model m- 

At maturity, the holder of a zero coupon bond receives one unit of currency. 
In accordance with the no-arbitrage principle, the time-f value of the zero-coupon 
bond with maturity 0 < t < T is modeled by 

P(t,T) := E(e~^ rhdh \T t )- (36) 

Translating this conditional expectation formally into an evolution problem of form 
©, we obtain g{x) = 1 as initial condition. We now have to realize that there 
is no weight 77 £ R d such that x K► £ L 2 ( R d ). We therefore split the 

initial condition into summands that each lie in a weighted L 2 -space. In the one¬ 
dimensional case, for example, we have g = l(_oo,o] + 1 ( 0 ,oop where l(-oo,o] G 
L 2 _(R) for every 77 “ > 0 and l(o,oo) £ L^+(R) for every 77 + < 0. 

Remark 5.2. We split the initial condition g into in 2 d summands gi that are 
supported in the 2 d orthants. To be precise, for j = 1,..., 2 d , let pi := (p{,... ,p J d ) 
with pj £ {— 1 , 1 } for the 2 d different possible configurations and let 

C)i := {(xi,..., Xd) £ R d | p{xi > 0 for all i = 1 ,..., d}. (37) 

By linearity of expectation, respectively of the PIDE, the problem can be split ad- 
ditively in 2 d separate problems. If for each of the summands gi a weight pi £ R d 
exists such that gi £ L 2 3 (R d ), then the results of Theorem \3.4\ can be applied to 
each problem with initial condition g J separately. 

As in Remark f5.2l we split the unity in the following way: 1 = g(x) = X}j=i lo-> ( x ) 
a.e. with the distinct orthants O 3 of R d given by m- For each j , we choose 

pi := —edr x i 2 pd (38) 

so that 1 pi £ L 2 ( R d ). If the distribution of Lt has a Lebesgue density, we 

may rewrite equation (ffiD as 
2 d 

u(T — t,x) = ^ u 3 (T — t,x) with u j (T — t,x) := E x ( 1 0 3 (Lt) e~ & rhdh y ( 39 ) 
j =1 


18 


K. GLAU 


Corollary 5.3. For e > 0 and a £ (0,2], let L be a time-inhomogeneous Levy 
process such that Ee e \ Lt \ < oo for every t < T and its symbol A satisfies (A2)~ 
(Af) for index a and every weight g £ R d with \g\ < e. Then for every 0 < t < T, 
the price of the zero coupon bond in model (|55|) is given as 

2 d 

P{t,T) = '$2u j (T-t,L t ) a.s., 

j =i 

where u 3 is the unique weak solution in W 1 (0, T; U“/“(R d ), LA (R d )) of 

v? + Ar-tu 3 + ru 3 = 0 , u(0) = 1 qj ■ (40) 

Proof. The assumptions yield that for each j = 1 ,..., 2 d , Conditions (A1)-(A4) are 
satisfied for weight rf and index a. According to Remark 14. 11 the distribution of 
Lt has a Lebesgue density, which yields equation (1551) . Now, the assertion follows 
directly from Theorem l3.4l □ □ 

It is worth mentioning that with the same technique we can characterize prices 
of options on a zero-coupon bond by solutions of PIDEs. A distinctive feature of 
the resulting PIDE is that the solution u of equation (1551) appears as the initial 
condition. Its initial condition thus is given by the solutions to PIDEs (15P1) . 

Interesting related applications are bankruptcy probabilities in the model of (Q), 
the value of barrier strategies in the bankruptcy model of (@) and reduced form 
modelling of credit risk as in (33). 

5.3. Laplace transform of occupation times. We characterize Laplace trans¬ 
forms of occupation times of time-inhonrogeneous Levy processes via weak solutions 
of PIDEs. Setting n := for D C R d Borel measurable, / := 0, initial condition 
g := 1, and inserting L 0 = x, equation (15U1) from Theorem 13.41 becomes 

u(T,x) = E x (e-' r -r° lD(Lh)d h), (41) 

which is the Laplace transform at 7 of the occupation time / Q T 1 o{Lh) d h that the 
process L spends in the domain D until time T. Landriault, Renaud and Zhou 
(2011) analyse Laplace transforms of occupation times of spectrally negative Levy 
processes using fluctuation identities. We characterize these transforms for a wide 
class of time-inhomogeneous Levy processes by parabolic PIDEs. Let us point out 
that the assertion is not restricted to spectrally negative processes as the examples 
of NIG and tempered stable processes show, see Section [4] 

Splitting the corresponding initial condition according to Remark 15.21 we let 

u 3 (T, x) :=E x (l oj (LT)e-^° lD ^ dh ). (42) 

Arguing as in the proof of Corollary 15.31 and applying Corollary 13.51 we obtain: 

Corollary 5.4. For e > 0 and a £ (0,2], let L be a time-inhomogeneous Levy 
process such that Ee e ^ Lt ^ < 00 for every t < T and its symbol A satisfies (A2)~ 
(Af) for index a and every weight rj £ R d with \r/\ < e. Let if := —ed~ 1 ^ 2 p 3 
as in (1551) . Then, u 3 from equation m is the unique weak solution in the space 
W 1 (0, T; if “/ 2 (R d ), LA (R d )) of 

u 3 + Ar-tu 3 + 1 du 3 = 0 , it( 0 ) = 1 0 j 

and u from equation (ED is given by 

2 d 

u(T, x) = ^2 ,u 3 (T,x). 

. 7=1 


(43) 
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If d = 1 and a > 1, then x i-» u(t,x ) := E x (e 7 ^o ^o(L h )Ah^ ^ er continu¬ 

ous with A = for each t £ [0,T] and in particular also continuous. 

5.4. Penalization of the domain. Observe that the limit of e _7 ^> 
as 7 —>• oo links occupation times to exit times. This idea lies at the basis of the 
repeated use of occupation times for modeling. Moreover, it opens a way to establish 
a Feynman-Kac-type representation of type (1201) for boundary value problems. In 
the language of diffusions, the presence of particles in the outer domain is penalized 
stronger and stronger until it is finally killed the moment it leaves the domain. For 
jump diffusion processes, the argument is outlined in | 6 ). In © and a forthcoming 
article, d2(i i, a similar technique is used for time-inhomogeneous Levy processes. 
Interesting for finance, the resulting Feynman-Kac-type representation serves to 
characterize prices of barrier and lookback options in pure jump models. The 
argument is based on the following result and the convergence of solutions for a 
sequence of killing rates of indicator type given by n x (x) := Al-Q«(a;) for A —» oo. 

Corollary 5.5. For a £ (0,2] and 77 £ R , let L be a time-inhomogeneous Levy 
process satisfying assumptions (Al)-(Af). Let f £ L 2 (0,T; if* (R d )) for some l > 0 
with l > {d — a)/2, g £ L 2 (R d ), k : [0,T] x R d —> R measurable and bounded, A > 0 
and D C R d open. Then the unique weak solution u x £ W 1 (0, T; H^ 2 (K d ), L 2 (R d )) 

of 

d t u x + Ar-tu x + KT~t.u x + Al-jjc M A = /, it A (0 )=g, (44) 

has for every t £ (0,T] almost surely the stochastic representation 
u x (T - t , L t ) = E(g{L T ) e“ £ Kh(ih -) dh e~ x £ 



Proof. The assertion follows directly from Theorem 13.41 □ □ 


5.5. Relativistic Schrodinger equation. Our analysis leads us back to the ori¬ 
gin of Feynman and Kac’s deep link between Schrodinger’s equation and diffusion 
processes. Recast in a relativistic mold, the formalism brings Normal Inverse Gauss¬ 
ian Levy processes into the spotlight: We find that a specific NIG process plays the 
same role for the relativistic Schrodinger equation as the Brownian motion does for 
the classical Schrodinger equation. Carmona, Masters and Simon (1990) provide 
a Feynman-Kac-type formulation of this link but give no formal proof. Baeumer, 
Meerschaert and Naber (2010) exploit this relation to model relativistic particle 
diffusion by an NIG process. We follow their presentation of the connection be¬ 
tween the relativistic Schrodinger equation and NIG processes. Then, Theorem 13.41 
allows us to make this link formally precise. 

The nonrelativistic Schrodinger equation for a single particle in a quantum sys¬ 
tem described by the potential energy V : R d x R + —> R is the following partial 
differential equation for the wave-function if : R d x R + — > C, 


= ( - ^ A + V ( x ^t)^ip(x,t), (46) 

where i is the imaginary unit, ^ denotes the time derivative of if, 2irh is Planck’s 
constant, m is the particle’s mass, and the Laplace operator A is given by Aif(x, t) := 

T.Ui f (*.<>• 

For a free particle, i.e. if V = 0, a formal connection to the Kolmogorov back¬ 
ward equation of the Brownian motion is obtained by the analytic continuation of 



20 


K. GLAU 


the Schrodinger equation fifjl) in time and inserting r = it. For V ^ 0, setting 
V(x,it ) := V(x,t) for every x and t, this relates equation (HH1) to 

h^(x,t) = (^~A -V(x,t)^(x,t), (47) 

which is the Kolmogorov backward equation of the killed Brownian motion with 
volatility a = \JH/2m and killing rate V/H. 

Let us now pass to the relativistic Schrodinger equation. According to Baeumer, 
Meerschaert and Naber (2010), the relativistic kinetic energy of a particle with rest 
mass m and momentum p is given by 

E(p) = i/||p|| 2 c 2 + m 2 c 4 — me 2 , (48) 

where c denotes the speed of light. The relativistic energy C%1) serves as a pseudo 
differential operator to define the relativistic Schrodinger operator 

«oW0M) (49) 

for the free particle. Thus, the relativistic Schrodinger equation for a single particle 
in a quantum system described by the potential energy V is given by 

i h^(x,t) = (Ho + V(x,t))ip(x,t). (50) 

Analogous to the nonrelativistic case, formally inserting r = it in equation m 
and setting V(x,it ) := V(x,t) for every x and t. yields 

^OM) + \(Uo + V(x,t))ip(x,t) = 0. (51) 

We note that j-E(p) is the symbol of the NIG process L with parameters a = me 2 , 
0 = 0,6=^, p = 0 and A = c 2 Id^, where we use the notation of Example 14.61 and 
Idd denotes the identity matrix in R d x R d . 

The following corollary formally establishes the Feynman-Kac-type relation of 
equation (l5ll) to NIG processes in terms of weak solutions. Note that here the 
potential energy V is allowed to be discontinuous. 


Corollary 5.6. Let the potential energy V be measurable and bounded. Let g £ 
L 2 ( R d ) for some rj £ R d such that \\ri\\ 2 < m 2 c 2 . Then the unique weak solution 

u £ W 1 (0,T- H^ /2 (R d ), L 2 (R d )) of 

u + -(Hqu + Vu) = 0, u(0)=5, (52) 

h 

has for every t G (0,T] the stochastic representation 


i(T-t,L t ) = E(g(L T )i 


-U t T Vt- h (L h )dh 


T t 


Proof. Corollary 15.61 is a direct consequence of Theorem 13.41 and Example 14.61 


(53) 

□ 

□ 


6. Numerical implementation 

Let us now explore the practical benefits of our Feynman-Kac-type result. We 
therefore implement a numerical scheme to solve Kolmogorov equation 0 for pric¬ 
ing path dependent options in jump models. We specify a class of employee options 
so as to shed light on the effect of a discontinuous killing rate. In order to give in¬ 
sight in the technique, as numerical scheme we choose the wavelet Galerkin method 
as developed by (13711 for European option pricing. This is a very powerful method, 
which uses compression techniques and can be adapted to more involved pricing 
problems, as we will demonstrate by incorporating killing rates. The implementa¬ 
tion requires some results from the classical theory on numerical analysis on partial 
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differential equations. In addition, the jump part of the operator needs some special 
treatment. We take care of the derivation of the discrete scheme by presenting the 
discretization steps (l)-( 6 ) below. 

We specify a type of employee option as described in Section 15.11 So as to 
include penalizations of low stock values permanently rather than only at a fixed 
maturity, we combine a call option with an indicator type killing rate. We specify 
the latter as instantaneous penalization for stock values below a fixed level by 
setting k(S) := —Al(_oo b](S) with a scale factor A > 0 and level B in equation 
ED. be. the payout at maturity is given by 

G(S T )e~ /o TA 1 <-~.B]( s *) dh ) 

where G(S) := max {5 — K, 0}. As driving process L in the model S = So e L we 
choose a pure jump Levy process from the family of CGMY processes described in 
Example 14.91 with parameters C > 0, G > 1, M > 0, Y £ [1,2), and whose drift b 
is given by the no-arbitrage condition (EFU) . Then, according to Example 14.91 the 
process and its symbol satisfy Conditions (Al)-(A4), with weight rj £ (G, —1) and 
index a = Y. 

We now fix a weight r) £ (G, — 1), denote by A the Kolmogorov operator of the 
process and let G(x) := G(e x ) and k(x) := —n(e x ). According to Corollary 15.11 
we obtain the fair price of the option by computing the unique weak solution u £ 
W 1 (0,T;^ / 2 (1R,),L 2 (R)) of 

u + Au + (r + k)u = 0, it(0) = G. (54) 

In order to prepare the discretization with finite elements, we first modify and 
then localize the equation to a bounded interval. The variational formulation of the 
resulting equation then allows us to discretize the space with a Galerkin method. 
Finally, the time discretization completes the fully discrete scheme. In more detail, 
we proceed along the following steps: 

(1) Modification of the equation: 

Choose a function if £ W 1 (0, T ; H^ 2 (\ R), L 2 (R)) such that cf := (u — tf) £ 
W 1 (0,T;IL i / 2 (]R.),L 2 (R)) and a;)| —»■ 0 for |ar| —»■ oo. Then (f is the 
unique weak solution of the modified equation 

<j> + A(f> + (r + k)(f> = /, 4>(0) = G - if(0). (55) 

(2) Truncation to a bounded domain: 

We localize the equation to a bounded interval (i?i,I? 2 ) with zero con¬ 
straints outside of the interval. Here, we for the first time encounter a 
conceptual difference between jump and non-jump processes: The jump 
part of the process renders the operator A nonlocal. It does therefore not 
suffice to specify zero boundary conditions. Rather, the values have to be 
set on the whole outer domain R\(I?i, I? 2 )- Formally, we incorporate these 
zero constraints by defining the solution space as ^ 2 {Ri 1 R 2 ) := {u £ 
H y G(R)\u\ [RuR2] c = 0} and L 2 (Ri,R 2 ) := {g £ L 2 (R)| 5 | [fll ^ 2 p = 0}. 
To be precise, instead of solving equation (l55l) we approximate the unique 
weak solution </> £ W 1 (0, T; U 1 ^ 2 (Ri, R 2 ), L 2 (R±, R 2 )) of 

4> + Af>+{r + k)4>= /, 4>(0) = (G-'0(O))l( flli _ R2) . (56) 

We now have to realize that we have changed the problem and that we 
need to control the resulting error \\<f> — f>\\ with an appropriate norm || ■ ||. 
Put differently, we have to choose the function if in step (1) in such a way 
that the error ||^ — <^|| decays fast as ~Ri,R 2 —> 00 . 
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(3) Variational formulation of the equation: 

Weak solution (j) G W 1 (0, T; iJ F / 2 (i?i, R 2 ), L 2 (i?i, R 2 )) solves operator 
equation dSHl) if and only if <j> satisfies the initial condition of (15B1) as a 
limit in L , that is 

lim^(t) = (G-V(O) t (RuR2) ) in L 2 (R 1 ,R 2 ) 

and for all v G Cq°(0,T) and ip G H y / 2 (Ri, R 2 ), 

T T 

~ J (4>(t), i P)L2i'(t) dt + j a(<j){t),ip)v(t)&t = f{<p,v), ( 57 ) 

0 0 

with bilinear form a : Ft' 1 / 2 (i?i, R 2 ) x H y / 2 (Ri 1 R 2 ) —> 1R, and f(ip,v) := 
fo f)l 2 v(t)dt. For the simplicity of presentation, we assume from 

now on that if is constant in time. 

(4) Space discretization with a Galerkin method: 

Coming to the heart of the Galerkin method, we choose a countable Riesz 
basis {w\,w 2 ,...} of H y ^ 2 (Ri, R 2 ) and define 

X n := spanjzci,..., w n } for all n G IN. 

Since H y / 2 (Ri, R 2 ) is dense in L 2 (Ri,R 2 ), we may further choose h n in X n 
such that h n —>■ <f >0 in L 2 (R±, R 2 ). We obtain the Galerkin equations for 
each fixed n G IN simply by restricting the variational equation (1571) . The 
resulting problem is: Find a function v n G W 1 (0, T; X n ; L 2 (i?i, R 2 )nX n ) 
that satisfies for all \ G Cq°(0, T) and ip G X n , 

T T 

- <p)l 2 x(t) dt + J a ( Vn (*)> v) x ( f ) dt = /(^> x) 

0 0 

V n (0) = hn- 


Elegantly, the classical theory guarantees the convergence of the sequence v n 
to <j> already in this abstract set ting . For more details we refer to Theo¬ 
rem 23.A. and Remark 23.25 in (fig ). 

The actual performance of the scheme, though, critically depends on the 
choice of the Riesz basis, which determines the rate of convergence. 

(5) Matrix formulation of equation (1581) .' 

Thanks to the linearity of the operators, we can simplify equation (1581) . 
Namely, it is enough to insert the basis functions w \,..., w n as test func¬ 
tions ip G X n in equation (l58l) . Then, denoting h n := '}2k=i a k w k and 
v n (t) := Y^ik=i Vk(t)wk, equation (l58l) turns out to be equivalent to 

n n 

T Vk(t)(w k ,Wj ) L 2 +'£Vk{t)a(w k ,w J ) = -a(ip,Wj ) 
k= 1 k— 1 

Vfc (0) = ctk for all & = 1 ,..., n. 


Written in matrix form the problem is to find V : [0,T] —» JR™ such that 

MV(t) + AV(t) = F (59) 

V(0) = a, (60) 
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where F = (Fi,...,F n ) T with Fk = —aytp,Wk) for k = 1 a = 

(«i,..., a n ) T , and the mass matrix M and stiffness matrix A are given by 

Mj k = (w k ,Wj) L 2 , Ajk = a(wk,Wj) for all j, k = 1 ,... ,n (61) 


Let us mention two critical points that arise in our setting. First, ap¬ 
proximation errors in the entries of the stiffness matrix A typically lead 
to significant numerical errors of the resulting scheme. As a consequence, 
they have to be computed with high precision. Second, due to the nonlo¬ 
cal nature of operator A , the matrix A is fully populated. This leads to 
a high computational cost of the solution scheme, which can be reduced 
considerably by using compression techniques. 


( 6 ) Time discretization: 

Having reached equations (l5f)l) and (16U1) , we are finally left to solve a linear 
system of ordinary differential equations. A variety of discretization meth¬ 
ods for these types of equations is available, for instance Euler schemes. 

To illustrate the numerical effect of the killing rate we use an implementation 
subject to the following specifications in the steps ( 1 )- ( 6)0 

• The equation is modified according to the choice ip(t, x) := maxje 1 — K , 0). 
Dominated convergence yields \u(t,x) — if>(x)\ —> 0 for |x| —>• oo. For a 
similar situation Proposition 4.1 in (1131 ) shows exponential convergence. 
We conjecture that also in our setting we have an exponential decay of the 
difference \<f> — <j>\. 

• As Riesz basis a wavelet basis of first polynomial order is chosen and com¬ 
bined with a compression technique replacing the stiffness matrix by a 
sparse one. We refer to Section 12.2.2 in the monograph of Hilber, Reich- 
mann, Schwab and Winter (2013) for a presentation of the wavelet com¬ 
pression technique and to (|451 ) for a related error analysis. 


• As time discretization an /ip-discontinuous Galerkin method is chosen as 
the initial condition is not differentiable and a scheme selecting more time 
points at the beginning is advantageous. For details we refer to Section 12.3 


In our numerical experiments we consider different domains on which the killing 
rate is active. Each of these domains is specified by a parameter B according to 
k(x) = Al(_ 00 i i 0 g (s))(x). Figure Q] depicts our results for B = K = 100 and B = 70. 
In both cases the maturity (in years) is set to T = 1 and the strike to K = 100. 
The parameters of the process are set to C = 0.01560, G = 0.0767, M = 7.55 
and Y = 1.2996. Notice that the scales of the graphs on the left side are chosen 
differently from those on the right side. 

We see that the killing rate of indicator type displays an effect in all of the 
considered cases. While the difference between the call and the employee options 
peaks around the level So = B, the killing rate affects prices globally with fast 
decay on both sides. The effects are stronger for higher scale parameters A. This 
leads to a monotone order of the price curves, the higher the scale parameter A, the 
lower the price when all other parameters are kept equal. 


^The author gratefully acknowledges Christoph Schwab and his working group 
for letting her use their code, which implements the Galerkin method for pricing 
European call options in CGMY-models. 
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European Call Option with killing rate 



Differences to European Call 



European Call Option with killing rate 



Differences to European Call 



FIGURE 1. Effect of a killing rate of indicator type on 
a call option price in a pure jump Levy model driven 
by a CGMY process. Top: B — K — 100. Bottom: 
B = 70. Left: payout of the call option along with the 
prices of the call option and the employee options with 
A = 1,10,100. Right: differences between the call price 
and the prices of the employee options. 


7. Robustness of the weak solutions 

We provide a robustness result that shows that small perturbations of the data 
/ and g and, more critically, of the bilinear form a only have a small effect on 
the weak solution of Kolmogorov equation ©• The result is crucial for the limit 
procedure in our derivation of the Feynman-Kac-type representation in Theorem 

HOI 

Let X H X* be a Gelfand triplet. For t € [0,T] and each n € IN 
let A 7 respectively At be an operator with associated real-valued bilinear form aj 1 
respectively a t . We introduce the following set of conditions. 
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(Anl) There exists a constant C\ > 0 such that uniformly for all n £ IN, t £ [0,T] 
and u, v £ X, 

max { |a"(w, v) |, \a t {u,v)\} < C\ ||ti||x|M|x- (62) 

(An2) There exists constants C 2 , C 3 >0 such that uniformly for all n £ IN, t £ [0, T] 
and u £ A', 

min{a"(it, u), a t {u, u)} > C 2 \\u\\ 2 x - C 3 \\u\\ 2 H . (63) 

(An3) There exists a sequence of functionals F n : L 2 (0,T; H ) —»• R + such that for 
all n G IN and u, v G L 2 (0, T ; #), both F n (u) —>■ 0 for n —> oo and 


T 


o 


- a t )(u(f),u(t))| dt < F„(u)||u|| L 2 ( 0iT;ff ). 


(64) 


Let /", f £ 

L 2 (0, T; Lf) with /" —>• / ffli z (0,T; X*) and < 7 ", </ € H with g n -A <7 in H. Then 
the sequence of unique weak solutions u n £ TT 1 (0, T; Y, iL) of 


Lemma 7.1. Let operators A and A n for n £ IN satisfy (Anl) (An3) 


u n + = /", «"(0) = g n (65) 

converges strongly in L 2 (0,T', X) n C(0,T; H) to the unique weak solution u £ 
W\0,T-X,H) of 

ii + Atu = f, u(0) = g. ( 66 ) 


Proof. Fix some n £ IN and let zt", u £ IT 1 (0, T; A, H) be the unique weak solutions 
of equations (1^51) and dSUl) and let w n := u — u n . Substracting equation (1^51) from 
m and inserting w n as test function yields for every t G [0, T], 


t 


J (w n (s),w n (s)) ds + 
0 



w n (s)) d s 


t t 

= J (f n (s)-f(s),w n (s)) ds + J (a" — a s )(u(s), w n (s)) ds. 
0 0 


(67) 


We insert f* (w n (s), in"(s)) ds = \{\\w n (t)\\ 2 H — ||u;"(0)||| f ), see e.g. |47t ) (equation 
(2) on p. 394), inequalities (1M1) . (1M1) and the inequality of Young. Subsequently 
applying the lemma of Gronwall yields the existence of constants c \, c 2 > 0 such 
that 

sup \\w n {t)\\ 2 H + Ci\\w n \\ 2 L ^ QT . x) 

t£[ °’ T] ( 68 ) 

< c 2 (\F n (u)\ + ||/" — /||i 2 ( 0) + \\9 n ~~ fflln) 

with F n from condition (An3). Hence u n u converges strongly in L 2 (0,T;AT) 
and in C(0, T; H), which proves the lemma. □ □ 


8. Proof of the Feynman-Kac-type formula, part (ii) of Theorem 13.41 

The key steps in the proof of the Feynman-Kac-type formula in Theorem 13.41 
are first applying Ito’s formula with the help of the regularity assertion in Lemma 
18.11 below, second invoking the convergence of regularized solutions to the solu¬ 
tion of Kolmogorov equation (O due to robustness result Lemma rm and third 
linking convergence in L 2 (lR d ) respectively L 2 (0, T; H^(R d )) to the convergence of 
conditional expectations via Lemma 18.21 
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Lemma 8 . 1 . For 77 £ R d and a > 0, let A be a pseudo differential operator 
whose symbol A has Sobolev index a uniformly in [0,T] x R v and let the mapping 
1 1 -)- At^—irf) be continuous for every £ £ R d . For k £ L°°([0, T] xR d ), g £ L 2 (R d ) 

and f £ L 2 (0,T;ff“ / 2 (R d )), let u £ W 1 (0, T; ff“ / 2 (R d ), L 2 (R d )) be the unique 
weak solution of 

u + A T-t.u + n T -tu = /, ( 69 ) 

«(0) = 9- (70) 

Then the following assertions hold. 

(i) Let m > 1. If g £ R d ), / e L 2 (0,T;^ m_ 1 )a/ 2 (R d )) and nh £ 

L 2 (0, T ; H^ a ' 2 (R d )) for all 1 < k < m and h £ L 2 (0,T-H^ a/2 (R d )), then 
u £ L 2 (0,T;i7™“ / 2 (R d )) and u £ L 2 (0, T; ^ m_ 2 )a/ 2 (R d )). 

( ii ) If g £ Hfj(R d ) for fj = m + d/2 + max(a, 1/2), / £ L 2 (0, T; H^(R d )) for 
7 = m + (d + l)/2 and k £ Co°([0,T] x R d ), then for every multiindex k = 
(fci,..., kf) with |fc| < m the derivative (1 + dt)D k u is in C([0,T] x R d ). If 
moreover, A is the Kolmogorov operator of a Levy process and f is continuous, 
then equation (17>§1) holds pointwise for all (t,x) £ (0,T] x R d . 


Proof. We derive the regularity assertion by explicit operations on the Fourier trans¬ 
form of the unique weak solution u £ W 1 (0,T; H^ 2 (R d ), L 2 (R d )) of equations (l69l) 
and CUD- We show the identity 

u = u := u 1 + u 2 + u 3 (71) 


with 

T,{u\t)) : .F,(„)e 

t 

T v (u 2 (t)) := I F v (f(s)) e-£ ,A *(-^) d »d s , 

0 

t 

J- V (u 3 (t )) := — j J,(fiM(s))e _ ^ T -‘ a a(--« 7 ) dA jg 
0 

and hence 

dt^niu 1 ^)) = -A T - t {■ - 

d t T v (u 2 {t)) = —A T - t {- - ir])T v (u 2 {t)) +.F„(/(£)), 
d t fF v (u 3 (t)) = -A T - t {■ - iri)F v (u 3 (t)) -T^nuft)). 


In particular, u satisfies equation (USD- Inequality ED from Remark mu with 
constants C± ,C 2 > 0 and the inequality of Cauchy-Schwarz guarantee the existence 
of constants C\, C 2 > 0 which are such that for all [t, £) £ [0, T] x R d , 


|-^(« 1 W)( 0 | <C'i|J-,( 5 )( 0 |e- tC2| « l “, 


1/2 


1/2 


\T v (u j (t))(0\ < c\ | J |^(/(s))(^)| 2 d S j e-^ 20 ^ 0 ds 

<c 2 (j |^(/ J '(s))( 0 | 2 (l + |^|) _a d S 


1/2 
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as well as 

|j' I) (5 t ii 1 (t))(0| < ci|^( fl )( 0 |(l + |£|)“e _tC ' 2|4r , 

1/2 

+ |W(*))(0|> 

for j = 1, 2 with f 1 = f and f 2 = —km. Hence there is a constant C 4 > 0 with 

L 2 (0,T-,H™ a/2 (R d )) + L 2 (0,T;Hi 1 m ~ 2 ' ,a/2 (R d )) 

— 1 )“/ 2 (ia<J))+ Wf\\L 2 (0,T-Hi l m ~ 1)o ‘ /2 (K d ))^~ L 2 (0,T-Hi, m ~ 1 ' ,a/2 (R d ))) ' 

For to = 1, by inserting u G L 2 (0, T; 7L)^ 2 (R d )) and km G L 2 (0, T; /f^ 2 (R d )), 
we obtain u G L 2 (0,T; H^ 2 (R d )) and d t u G L 2 (0, T; ff/~ a,/ 2 (R d )). In particular, 
u G W 1 (0,r ; ^ /a (R d ),L2(R d )) is the unique weak solution u = u of equations 

AMD and (1701) . 

For to = 2 it is thus sufficient to notice that km G L 2 (0,T;IL^ 2 (R d )) implies 
m G L 2 (0,T; H^ 2 (R d )) and d t u G L 2 (0, T; L 2 (R d )). An iterative argument then 
yields part (i) of the Lemma. 

(ii) By the inequality of Cauchy-Schwarz and f Rd (l + |£|) d e d£ < 00 if e > 0, 
we obtain that for /3 = to + d/2 + max(a, 1 / 2 ) and 7 = to + (d + l)/2 there exists 
a constant C 5 > 0 such that 

/1(1+ ^((m 1 +u 2 )(t)) (oi (1+ieir ^ 

R d 

< C 5(llslli^(Ry) + ll/lli 2 (0,T;^(R d ))) < °°- 

Furthermore, the mappings t >-»• (u(t)) (0 and t i-» (&(£)) (£) are con¬ 

tinuous for each £ G R d . Dominated convergence implies D*( 1 + 9 j)(m 1 + m 2 ) G 
C([0, T] x R d ) for every multiindex k = (Aq,..., fed) with |fc| > 0. Moreover, there 
exists a constant Cq > 0 such that 

f |(l + <9 t )j^(u 3 (f))(£)|(l + |£|) m d£ ^cellKtill^o R d )) < °°' 

R d 

Dominated convergence yields 1 + d t )u 3 G C([0, T] x R d ) for every multiindex 
k = (fci,..., kd) with |fc| > 0 . 

Now let A be the Kolmogorov operator of a time-inhomogeneous Levy process. 
In order to establish equation (l69l) pointwise, fix a t £ T for which the equation 
holds (as operator equation) and choose a sequence u n G Cq°((0,T) x R d ) such 
that u n (t) —> u(t ) in the norm of H^ 2 (R d ) and u n (t) —> it in the norm of L 2 (R d ). 
Moreover, let <p G C'“(R d ). We notice that Ar-tu(t) is defined pointwise, since 
u(t) G C 2 (R d ). An elementary manipulation and the continuity of the scalar prod¬ 
uct yield 

( A T -tu{t 1 x)ip{x)e~ 2( ' ri ’ x) &x = (u(t),A^*<p ) L 2 = lim (m„(4^/!)V)i 2 

I 77 n—>-oo 77 

R d 

with the adjoint operator A ///* defined in Lemma lA.2l in Appendix 1X1 Equation ([7]) 
from the same lemma and the continuity of the bilinear form imply 

lim {u n (t),AT-?<p ) L 2 = lim a T -t(u n (t),ip) = a T -t(u(t),ip) 

n—>• 00 77 n—>• 00 


|j- ?? (a t M j (i))(^)| < c 2 
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and hence 

(u(t),ip ) L 2 + {A T -tu(t),<p ) L 2 = (f(t),ip ) L 2 for all <p £ (!R d ). 

From the fundamental lemma of variational calculus, u(t, x)+Ar-tu(t, x ) = f(t, x) 
follows for a.e. x £ R d . Since we can choose t arbitrarily from a dense subset in 
( 0 ,T), the assertion follows by continuity of u + At—U — /• □ □ 


Lemma 8.2. For rj £ R d and a £ (0, 2], let L be a time-inhomogeneous Levy pro¬ 
cess with symbol A = (^4t)te[o,T] that satisfies exponential moment condition (Al) 
and Carding condition (A3). Then 

(i) for every t > 0 there exists a constant C(t) > 0 such that .E[|</j(L s )|] < 
C'(t)||<^||i 2 (] R d) uniformly for all ip £ L 2 (R d ) and s £ [t,T\, 

(ii) for l > (d — a)/2 and every 0 < t < T there exists a constant C\ > 0 
such that |i?( f f T <p(s, L s ) ds| iFt) | < Ci\\tp\\ L 2 ( t T . H i ( R d)) uniformly for all tp £ 
L 2 (0,T;^(R d )). 

Proof, (i) By Remark H~Tl and Condition (A3), the distribution of L t has a Lebesgue 
density. Applying Parseval’s identity, we obtain 

E H L t)\ = J T(\ip\)(f~iri)e-ti A ^-^ d sdf. 

n d 

Inserting inequality m and the inequality of Cauchy-Schwarz then yields asser¬ 
tion (i). 

(ii) W.l.o.g. ip > 0. We have E( f { T <p(s, L s ) ds| = G(L t ) with 

T-t 

G(y)=E^ J tp(s + t, L t + S - L t + y) ds). 
o 


The theorem of Fubini and Parseval’s identity imply 

T-t 

G(y) = J J tF(r y ip(s + t)) (£ - irj)e~ fo A t+ ^-ir,) d« ds 

R d 0 


where T y f(x) := f(x + y). Notice that T ri (r y f )(£) = e T , ? (/)(£). Inserting 

the inequality of Cauchy-Schwarz and equation m with constants Ci, C 2 > 0, we 
obtain that for l > d — a there are constants ci, C 2 > 0 such that 


' T-t 


T-t 


1/2 


\G(y)\ < Ci J | J \T(T y ip(s + t))(^-i v )l ds J e 2sC2| « |Q dsj d£ 


R d \ 0 
T 

<ci// \iF(<p(s + t)) (£ — ip) | (l + 

0 R d 

— C2 ^P^L 2 (0,T-H l f 2 (R d ))' 


d£ ds 


This concludes the proof. 


□ 


□ 


We are now in a position to prove part (ii) of Theorem 13.41 

Theorem \3.4\ part (ii). First, assume that t n- At(£ — ip) is continuous for all 
£ £ R d . By density arguments, respectively mollification, we can choose sequences 






FEYNMAN-KAC FORMULA FOR LEVY PROCESSES 


29 


g n G C£°(K d ), /” e C£°([0,T] x K rf ) and re” G C£°([0,T] x K rf ) such that for 
n —>• oo, 

g n ^g inLl(R d ), 

f n ^f inL 2 (0,T;^(R d )), 

re” —>• re pointwise and sup ||re”|E°o([o,T]xR' i ) < oo. 

raeW 

We denote, with a slight abuse of notation, by a” the bilinear form associated with 
At + re” and a* the bilinear form associated with At + Kt- Then 

at(u,v) = a t (u,v) + ((K™ — K t )u,v) L 2 ( R d' ) for all u, v G H^ 2 (R d ). (72) 

Together with Conditions (Al)-(A3) and the uniform boundedness of re” and re, we 
obtain the validity of Conditions (Anl) and (An2) from Section 6. Moreover, by 
equality (EH) and the Cauchy-Schwarz inequality, we get 


(a” — a) (u(s), v(s)) ds 


— ||( K K ) w || L 2 (0,T;L 2 (JR. d )) L 2 (0,T-,L^(Ii d )) ■ (73) 


It follows from pointwise convergence of re” —> re and dominated convergence that 
for n —» oo, 

F n (u) := ||(re” - «H| L 2 ( 0 )r . L 2 (R<i)) -t 0 for all u G i 2 (0, T; L 2 (R d )), 

and hence Condition (An3) is satisfied. Let it” G W 1 (0, T; H^ 2 (R d ), L 2 (R d )) be 
the unique weak solution of 


ii n + A t u n + n?u n = /”, u”(0) = g n . (74) 

Lemma 17711 yields the convergence it” —► it, both in the space L 2 (0, T; H^ 2 (R d )'j 
and in C(0, T; L 2 (R d )), to the weak solution u G W 1 (0, T; H^ 2 (R d ), L 2 (R d )) of 

u + Atu + Ktu = f , i/,(0) = < 7 . (75) 


Lemma 1701 shows that equality (I7H) holds pointwise and that it” is regular enough 
to apply ltd’s formula. We denote by (b t: <r t , F t \ h) t(z ^ Q „ the characteristics of L 
and set w n (t,x ) := u n (T — Lx). Then Ito’s formula for semimartingales, see for 
instance Theorem 1.4.57 in (13211 . entails 


w 


‘(T, L t ) e" 7o T “-I ( L 0 dA l s ) e“ % K Z( L ^ dA 


w n — A T w n — KW 71 


(r, L T )e~Jo K Z( L x)d\ dr 


1 

+ j (o-y 2 ■ Vic”(r,L r )) e-^^l^dlF, 

S 

+ (V7o K X (i *) dA (w”(.,L._ +x) - w”(-,L._))l (Sj00 )(-)) * (At - t/) T . 


(76) 


Thanks to our assumptions on g n ,f n and re”, we may decompose u” in three 
summands as in equation EH. Then, applying of part ( ii ) of Lemma 18.11 it is 
elementary to conclude that w n and Vic” belong to L 2 (R d ). Hence, the inte gral s 
with respect to W and g — v are martingales, compare Theorem II.1.33 a) in (32|). 

We insert the identity to" — A T w n — re”iu” = /" with f(t,x ) := /”(T — t,x) in 
equation ESI), subsequently multiply the equation with the term e^o K U L O d A anc j 
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lastly take the conditional expectation. This gives us for 0 < s < T, 
E (w n (T, L t ) e“ ft K " (ix) dA 

T 


F s ')-w n (s,Ls) 

= eU J n (T,L T )e-ft K ^ dX dT\E s ) . 


(77) 


Let w.l.o.g. 0 < s < T. We will now derive the desired stochastic representation by 
letting n —> oo for each term in equation (1771) . Denote w(t , x) := u{T — t , x). From 
the convergence w n (s , ■) —>• w(s, •) in L 2 (R d ) and part (i) of Lemma [FT?1 for s > 0, 
we get the convergence 

w n (s,L s ) —> w(s,L s ) in ld(P) and a.s. for a subsequence. 


The pointwise convergence and the uniform boundedness together with 

dominated convergence imply both f a K\{L\) dA —> f a K\(L\)d\ and uniform 
boundedness of the sequence for 0 < a < b < T. 

Together with w n (s,L s ) —>• w(s,L s ) in L 1 (P) the convergence 


E(\w n (t,L t )e 


— Sj K \(Lx)dX _ 


w(t, L t ) e 


-if «*(£*) d* 


E s 


as n —> oo then follows from the triangle inequality. 

Next, denote f(t,x) := f(T — t,x). Since /"—>•/€ L 2 (t,T; Lf^(R d )), part (ii) 
of Lemma (S21 guarantees the existence of a constant C 2 > 0 for l > (d — a )/2 such 
that 


T 

E( j \(T - 7 )(h,L H )\dh\T s ) 


II f n -FU 


Now the triangle inequality yields the convergence of the second line in equation (EB 
and therefore part (ii) of Theorem 13.41 under the additional assumption that the 
mapping t n- A t (£ — ir]) is continuous for every £ £ ]R d . Finally, thanks to the tower 
rule of conditional expectation, the claim follows by induction over the continuity 
periods also under the more general assumption that t At(£ — ir]) is cadlag for 
every £ G R d . □ □ 
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Appendix A. Adjoint Operator 

For a Levy process L with characteristics ( b , a, F: h) we denote by 
A ( b,a ’ F ) an( j j[(b,<r,F) p- s Kolmogorov operator and its symbol, respec¬ 
tively. As the following assertions straightforwardly extend to the time- 
inhomogeneous case, we here present the time-homogeneous case only. 
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Lemma A.l. For p G R d , let L be a Levy process with characteris¬ 
tics (b, cr, F\ h) that satisfies exponential moment condition EM{j]) and 
denote by A its Kolmogorov operator with symbol A. Then 

^(0 := A(f + ip) = "**>(€) + A(ip) for all £ G R d 

with 

b v — b + a ■ p + j (e {r, ’ y) -l)h(y)F(dy), 
n d 

Ffidy) =e {v ’ y) F(dy). 

In particular, A v is the symbol of a Levy process with killing rate A(ip). 
Moreover, its Kolmogorov operator A v satisfies 

A v p = e~ M A (e^ p) = A {bV ' a ’ FTI) p+A(ip)p for all ip G C™{R d ). 

Proof. It is elementary to verify the assertion for the symbol. This 
can then nicely be used to verify the assertion for the operator: Let 
<p G C' 0 oo ( R d ), then ^(e^ <p)(£) = E(p)(£ - ip) and 

A (e^ ip) (x) = J 7-1 (AF( e <^> p)) 

= Wfj e ’‘ K,I> “ ir >) d « 

R d 

= Wy> / + 

R d 

which concludes the proof. □ □ 

For all p G Co°(R d ) let 

Efip) := e-^-) E(pe^) an d E~\p) := e"^ E^(pe M ). 

Theorem 4.1 in (16) shows that for a pseudo differential operator A 
whose symbol A has a continuous extension to U- v that is analytic in 
the interior of U- v and satisfies the continuity condition (A2), we have 

A<p = J F-\Ap) = T~ 1 {A~ y Tfip)) for all <p G C 0 °°(R d ). (78) 

Parseval’s equality yields for all p, G Co°(R d ), 

a(p,fi) = (Ap ,^) L 2 = -^(A- Tl f ri {p),E r] {fi>)) L 2 . (79) 

For the pseudo differential operator A and its symbol A, we define 
their Lf ? -adjoints A v ’* and A~ v '* such that for all p, G (^(R^), 

(A,1p)L2 = (tp,A’ , -’4i) L ; t , 


( 80 ) 

( 81 ) 
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Lemma A.2. For r/ G R d , let L be a Levy process with characteris¬ 
tics (b, cr, F\ h) that satisfies exponential moment condition EM{rf) and 
denote by A its Kolmogorov operator with symbol A. Set F~ V, *(B) = 

f or Borel 3ets B 7 1 I 0 ). where F ;yV B ) = ^F-1(B) + 

F-'(-B) and FJ ym (B) = F-(B ) - F n l(B) and let &-"•* = -fr’. 
Then 

A~ v ’* = A^ = A + A(-iri), 

A ^V = A ( (p) = A ^^ + A(-irf)(p. 
Moreover, F _T? ’* is a Levy measure. 

Proof. For every ip G (7^(11^) we have 

(A-'iF y ( v ),F v (tj,)) L , = (J-V(eW <p),r{e^ tp))v 

= (F(e (ri -' } ip)) L i. 

Since A(z) G R for z G C d with 9ft(z) = 0, by Lemma [A. 1 1 we obtain 
~A~h = A( b ~ v ' a ’ F - r >) + A(-ir)). 


As A ( b ~ v ’ a ’ F ~”) is the symbol of a Levy process, we have 

A (*-W”)(f) = for all £ G R d , 

whence the assertion of the lemma follows directly. □ □ 


Appendix B. Proof of Theorem 13.31 

Theorem B.l. For r) G R d and a G (0,2], let the symbol A = 
(A t ) t& [ 0 ,T] of pseudo differential operator A = (A t )tG[o,T] be such that 
for each t G [0, T], A t has a continuous extension on U_ v that is ana- 

o 

lytic in the interior U_ v and there exist constants C(t),m(t) > 0 with 

\A t (z)\ <C(t)(l + \z\) m{t) for all z G U- v . (82) 

Then the following assertions are equivalent. 

(■ i ) The operator A is parabolic with respect to (//“, (R d ), L^, (R d )) 
uniformly in [0, T] x R v . 

(ii) The symbol A has Sobolev index 2 a uniformly in [0, T] x R rr 


The proof of the theorem is a straightforward generalization of the 
proof Theorem 3.1 in (2ll). where the assertion is proved for symbols 
that are constant in time and for spaces without weights, that is if 
ij = 0. In order to provide a self-contained presentation we give a 
detailed proof. 
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of Theorem \B.l[ By the assumption on the analyticity of A and in¬ 
equality (1821) . we obtain from Theorem 4.1 in (fl 6 h that for all t G [0, T], 
rf G and ip, if G C'^°(R d ), 

= J Mt- imnm - - irf) d£. (83) 

R d 

This equality yields that (Cont-A) implies (Cont-a). Together with 
the following elementary inequalities, it also yields that (Gard-A) im¬ 
plies (Gard-a): For C\ > 0, C 2 > 0, 0 < j3 < a and 0 < C 3 < C\ there 
exits a constant (7 4 > 0 such that C\X a — C 2 x^ > C 3 x a — (7 4 for all 
x > 0 and 

c 2 kl 2 " - C 3 (l + If I 2 )' 5 > C 2 |f I 2 " - C'(l + IfU) > c 2 (l + |f|) 2 “ - c 3 

(84) 

with a strictly positive constant c 2 and C 3 , c 3 > 0. 

Moreover, piecewise continuity of t ha a t (u, v ) for every u, v G Hf(R d ) 
follows from the piecewise continuity of t i —> A t (z ) for every z G U_ v 
and dominated convergence, which applies thanks to (Cont-A). 

For the implication from (i) to (ii), we first show the following. Let 
7 : H d —» R be a continuous function. If we have 

J 'y{Z)\Fr/(u)(Z)\ 2 e- 2 «’t ) dS>0 (85) 

R d 

for all u G i7“,(R d ) for which R r f(u) is compactly supported, then 
7(0 > 0 for all £ G R d . To prove this claim, we follow closely the 
derivation of the fundamental lemma of variational calculus. Let us 
for a moment assume that y(£) < 0 for some £ G R d . Due to conti¬ 
nuity, 7 is negative on a nonempty open subset of U C R d . We now 
choose a function u G H“,( R d ) such that its Fourier transform T v fu) is 
smooth, nonconstant and compactly supported in U. For this choice of 
u, however, the integral in inequality (1551) would be negative, leading 
to a contradiction. This shows that 7 > 0. 

We observe that (Cont-a) implies inequality (1551) for the continuous 
mappings £ ha (7(1 + |£|) 2q ± 3?(A f (£ — irj')) and £ ha (7(1 + |£|) 2 “ ± 
3(kb(£ — ir/)) for all t G [0,T] and rf G R v . Thus (Cont-A) follows. 
Similarly, using once again inequality (154)) . we obtain that (Gard-a) 
implies (Gard-A). 

Finally, we observe that lim^t a s (u, u) = a t (u,u ) implies that 

1 ini f 7f s (£-i?7 , )|7 r ?? '(u)(£)| 2 d£ = J A(£ - irf) (u)(£) | :2 d£, 

R d R d 
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while on the other hand dominated convergence shows that 
lim [ A s (£ = [ lim A s (f - irf) | 2 df. 

s—>t J J 

R d R d 

Now, an application of inequality (1851) yields lim,,^ A s (£—irj') = A t (f — 
ir)') for all £ G R d and rf G R r] ■ Therefore piecewise continuity of the 
bilinear form entails piecewise continuity of the symbol. □ □ 


Theorem 13.31 now follows as a Corollary from Theorem IB. 11 


Theorem [Of According to Lemma 2.1 (c) in (16|), for every 0 < t < T 
the map z K > A t (z) has a continuous extension to the domain 


Moreover, Theorem 25.17 in (1431) together with Lemma fA.il in Appen¬ 
dix El shows that inequality (182|) is satisfied for each t G [0,T] with 
m(t) = 2 and some constant C(t) > 0. Theorem 13.31 is now applicable 
and yields the corollary. □ □ 
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